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SUMMARY 


This investigation was performed at Princeton University in the 
period from February, 1951 to June, 1951. 

The original intention was an experimental investigation of the 
feasibility of using a powered model in harnessed circling flight to evalu- 
ate the longitudinal stability of a full-scale airplane. The airplane 
selected was the Cessna 140 two-place light airplane, since its stability 
derivatives and other parameters had been determined and were at hand. A 
one-tenth scale flying model was constructed. 

The method of attack was to have been as follows: 

A. Determination of stick-fixed neutral point, N, from steady- 
state flight tests. 

B. Experimental determination of the model's frequency response. 

§ by obtaining a transient response of © to an elevator pulse input 
and converting this to the frequency response by the method outlined in 
Reference ll. 

C. Correlation of A and B with theory. 

No quantitative data were obtained for A and B due to mechanical 
difficulties and lack of time. Sufficient qualitative data were obtained, 
however, to show agreement between the model's transient response and that 
predicted by theory; to point up several difficulties peculiar to this type 
of testing; and to indicate recommendations for further development of the 


method. 


It was found that the scale chosen gave too emall a model from com 


siderations of required mechanical installations and dynamic similarity. 
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Transient response tests of relatively short duration were found 


to be the only feasible type due to the presence of wake disturbances. 
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INTRODUCTION 


With the increasing number of new and radically different configura- 
tions of airplanes being introduced for sonic, transonic and convertaplane 
purposes it seems that some method of evaluating the low-speed longitudinal 
stability characteristics in unaccelerated forward flight is desirable prior 
to fixing a design. This could most easily be done in model form. A tech- 
nique such as this would allow the inertia forces on the model to be inter- 
preted to full scale results and would also be valuable in determining gust 
loadings. 

The Air Materiel Command at Wright Patterson AFB established a test 
facility which gave promise of yielding such data. 

An analytical study (Reference 5) was made by AMC. From this it 
was concluded that certain of the longitudinal stability characteristics 
could be evaluated, and that even where the dynamic stability of the phugoid 
mode was greatly modified by the harnessing forces on the model, the model 
response could be related to full scale response by certain mathematical 
operations. Unfortunately the test facility was decommissioned before any 
data were taken to substantiate the above conclusions, and at the time of 
commencement of this investigation, had not been re-established. 

This facility made use of powered models flown in harnessed circling 
flight; a type of flight known as "control-line flying" which is used as a 
sport or hobby by certain of the general public. 

It was the purpose, therefore, of the present investigation to 
provide a limited amount of information on which to evaluate harnessed cir- 
cling flight of models as a testing method, and to establish, if possible, 


testing techniques which will compleMent the method. To this end a one-tenth 
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scale model of the Cessna 140 light two-place airplane constructed. Reasons 
dictating choice of this airplane were that complete static flight tests and 
wind tunnel tests of the Cessna 140 have been previously made here by the 
Aeronautical Engineering Department of Princeton University and therefore 
complete data were available. This investigation was limited to comparing 
the neutral point, as obtained from level flight testing of the model, with 
the neutral point from flight test of the full-scale airplane; hall comparing 
the frequency response of the model's pitch angle to that obtained from the 
equations of motion of the airplane. The experimental frequency response 
was to be obtained from the analysis of the transient response to a triangu- 
lar pulse of elevator deflection. 

These limitations were dictated by lack of facilities, funds and 
time. 

The data were to be obtained from static stability flight tests of 
the model with three center-of-gravity positions and fram transient response 
tests at one center-of-gravity position (27.8% mac). 

It is believed that this investigation is the only one of its kind 
to date. Although the final results may not be conclusive due to rather 
primitive equipment, it may indicate a method and provide a part of the in- 
formation necessary for a complete analysis of the method, which of course, 


requires numerous tests on many more of the longitudinal stability aspects. 
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SYMBOLS AND ASSUMPTIONS 


AR - Wing aspect ratio 


AC - Aerodynamic center of the wing 








a & (ses) - Slope of the lift curve of 
AX 
: the wing (two dimensional) 

(per degree) 

aw = (<<) - Slope of the Lift curve of the 

Aha fw 

wing (three dimensional ) 
(per degree) 

Q, = (so) - Slope of the lift curve of 

Aa | + 


the tail (three dimensional) 


(per degree) 


b 


C - Mean Aerodynamic Chord of the Wing -- feet 


Wing span -- feet 


C.G.- Center of Gravity 


7 WPVv2s” . « Wire drag coefficient 
Cy = 2:+Dw - Drag Coefficient 

. “Pp ViS 
Co- Effective drag coefficient (C,-C-) 

an 

C=- - Thrust coefficient 

i 4P v*S / 
C, - Parasite drag coefficient 


C,- Induced drag coefficient 


©) = as 
E- as - Lift coefficient 
Cy, = = -~ Lift coefficient of the tail 
Cm: _ - Pitching Moment coefficient 
qSce 


of the airplane 
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Cy. = ACo - Change in drag coefficient 
Aan 

with change in angle of 
attack. 

C E AC M4 - Change in effective drag coefficient 

Con § 
AX 


with change in angle of attack. 


A > - Angle between flying wires and datum plane 








Change in lift coefficient 
with change in angle of 

attack -- per radian 

Change in pitching moment 
coefficient with change in 
angle of attack -- per radian 
Change in pitching moment 
coefficient with rate of change 


of angle of attack with respect 


Change in pitching moment co- 
efficient with rate of change 


of pitch angle with respect 


Change in pitching moment co- 


efficient with change in 


Cul = AC. - 
Ax 
- ae 
se 
Ge, a (oe = 
dx = 
d (= 
to ie 
Cinde = at : 
A= 
% 
# 
to In 
elevator angle. 
D - Total drag of the airplanes -- lbs 
d - Differential operator 
C - 0.869 - Airplane efficiency factor 
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med 
= 


F 
£ 
S| o 
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Centrifugal force in level flight 
‘2 
aay - Non-dimensional line force 
RpvVS 
Acceleration due to gravity -- 32.2 ft/sece 
xu 
A (ks) - Non-dimensional inertia parameter 
pec” 
: - Airplane moment of inertia about 
the Y axis -- slug ft 
Vol 


Mean aerodynamic chord -- per cent 

Radius of gyration about the Y axis -- ft 

Total lift of the airplane -- lbs 

Tail length (airplane CG to AC horizontal tail) -- feet 
Moment about the { axis -- lb-ft 

Mass of airplane -- slugs 

Miles per hour 

Stick-fixed neutral point 

Normal acceleration (along the Z axis at the CG) -- "G" wits 
Period of oscillation -- seconds 

Ze pr y* - Dynamic pressure 

Radius of flight circle. 

Wing area -- fte 

Horizontal tail area -- ft© 

Time to damp to 1/2 amplitude -- sec 

Time -- sec 

Speed ratio = » incremental change in airspeed divided 
by airspeed 


True airspeed -- ft/sec 


i 
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OS Se Lt 
V Ss ¢ 
\V - Gross weight (average) -- pounds 





- Volumetric coefficient of the tail 


GREEK SYMBOLS 

QM - Angle of attack of wing -- degrees 

A, - Angle of attack of the horizontal tail -- degrees ( 
dS, - Elevator deflection -- degrees 

€ - Angle of downwash -- degrees 

Np = Os - Airplane tail coefficient 


q 
0 - Angle of pitch (angle between % axis and horizontal -- degrees ) 


iS Root of stability quartic 
i - Relative airplane density 
proc 
factor (non-dimensional) 


fP - Mass density of air -- sluge/ft> 
7 
,T .-—-—_—_ - Time parameter -- seconds 
PSV 
() - Phase angle -- degrees 
GQ) - Angular frequency -- radians/sec 


G - Angular frequency, radians per t Vie 


A variable with subscript zero indicates the maximum value of that 
variable during an oscillation and may be either positive or negative. 


Vertical bars indicate the conventional absolute value of any variable. 


SIGN CONVENTION 
Reference Axes 
Normal stability axes (wind axes) are used. 
X-axis - longitudinal axis - positive forward 
Y-axis - lateral axis - positive along right wing 


Z-axis - vertical axis - positive down 
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Linear Displacements 


A linear displacement along a positive reference axis is considered 





positive. 


Angular Displacements 
An angular displacement which is clockwise when viewed from the 


origin looking along a positive reference axis is considered positive. 


Velocities and Accelerations 
Velocities and acceleration, either linear or angular, are con- 


sidered positive in the same sense as the corresponding displacements. 


Control Surface Deflections 
Positive elevator angle is associated with a downward movement of 


the elevator trailing edge. 


Phase Angles 
A negative phase angle signifies a time lag of the response to 


the forcing function. 
A positive phase angle signifies a time lead of the response to 


the forcing function. 


ASSUMPTIONS 
In the theoretical analysis and in the evaluation of the stability 
derivatives the following assumptions are made: 
1. Stability axes (wind axes) are used. 
2. The airplane is initially in level flight. (This means for 
stability axes cosH6=!1, sinQ@=4 ) 
3. All disturbances are small, so that the equations of motion may 


be considered linear. 








LO. 
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The lift caused by the horizontal tail is negligible. 
Small lateral disturbances which may exist do not affect 
the aerodynamic forces or pitching moments in the plane of 
symmetry. 
Forces due to pitching velocities are negligible. 
Change in pitching moment due to the rate of change of 
elevator angle is negligible. 
Only longitudinal motion is considered. The moment of inertia 
is that about the Y axis. 
Criteria for dynamic similarity were based on equal C, of 
model and airplane (References 5 and 10). 
Standard model flight conditions used in the solution of the 
equations of motion were: 

Y= 58.7 {pa 

R = 37-5 fk 
at datum level = 4 ft. above the floor. 


10 
181 
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DESCRIPTION OF THE MODEL AND INSTRUMENTATION 


The model used in this investigation was a one-tenth scale model 
of the Cessna 140 two-place light airplane. It was constructed by the 
authors. The construction consisted primarily of balsa-wood, paper-covered 
and doped; a construction similar to spin tumnel test models (Reference 9). 
Original intentions had been to make the model dynamically similar on the 
basis of equal lift coefficient (References 5 and 10). However » the choice 
of airplane and the scale used soon indicated that this would be physically 
impossible with available —— » due to the low weight requirement and 
small size of the model. The choice of airplane was fixed by considerations 


previously mentioned, and scale was fixed by (1) budgetary reasons plus 


(2) the requirement of an indoor flying site to reduce stray wind disturbances 


on the model. It was decided to stay as close as possible to the weight 
limitation and measure the resulting moment of inertia of the model, using 
that and the actual weight in the non-dimensional equations of motion of the 
airplane of the harnessed model. This modification in no way affects the 
ability to judge the effectiveness of this type of flight as a test method. 
The model was powered by a Bantam .19, cubic inch displacement, 
1/10 HP (@ 9000 RPM), internal combustion, air-cooled, model airplane engine 
with glow ignition (Figure | ). This ignition saved weight by requiring no 
ignition components in the model; merely a 1$-volt battery as a starting 
device. The engine was mounted inverted to preserve the contours of the 
upper cowling and it projected out of the lower cowling a short distance 


(Figure 2). A pressure type fuel feed system of commercial make was used 
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to prevent variation of engine speed due to change of fuel head in flight 
(Figure | ). | 

In order to give acess to the mechanism, engine and fuel system 
the wing, as well as various hatches were made removable (Figure $). 

In order to clarify the operation of the elevator pulsing mechan- 
ism it is necessary to describe the control of the model in flight. In 
normal flight the model flies in a circle around the operator on two wire 
lines which lead from a bellcrank (Figure 4) in the fuselage, through a guide 
: near the wing tip handle in the hand of the operator. The model is flown 
with some yaw to the outside of the circle which is produced by offsetting 
the rudder and by attaching the bellcrank behind the C.G. The amount of 
this yaw is also controlled by adjusting the line guide on the wing. A push 

elevator horn 
rod from the bellcrank links with the so that a motion of the operator tend- 
ing to increase tension in the top line from the handle and decrease tension 
in the bottom line is transmitted as an up deflection of the elevator. 

The requirements on the elevator for the pulse input are that it 
(1) deflects from the trim elevator position, (2) returns to this position, 
and (3) remains locked during the ensuing transient motion. To this end a 
third wire from the operator to the model was employed. When a level equili- 
brium flight condition was attained, the operator transferred tension from 
the two flying wires to the third wire. The point of attachment of this 
wire was slightly aft of the model's C.G. in order to maintein slight outward 
yaw and hence line tension. The first result of tension on the third line 
was the application of a brake on the elevator push rod (Figure 4). ‘This 


locked the elevator in the equilibrium position previously established. As 
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the tension was increased to counteract the full centrifugal force of the 
model, the third line also triggered a rubber-powered pulsing mechanism 
(Figure 4 ). The cam in this mechanism displaced the aft portion of the } 
push rod, working against a spring which returned the elevator to the original 

locked position. The cam was cut to give a pulse of triangular shape with 


time, having a duration of approximately .04 seconds. By adjusting the posi- 





tion of the cam follower on the push rod, the amplitude of the pulse was 
varied. The cam contacting the follower closed an electrical circuit through 
a third wire which fired a flash bulb. After the elevator had been pulsed, 
the mechanism was automatically re-cocked for another rm when the tension 
on the third line was released and the elevator again controlled by the oper- 
ator. 

Measurement of elevator angle was necessary in both phases of the 
testing. In the steady-state tests the angle required for equilibrium at 
each velocity had to be known, and in the dynamic tests the determination of 
the exact size of the pulse was necessary. This angle was to be measured 
with a Giannini Micro-torque. This instrument is an electrical potentiometer 
employing a rotating contact of very low mechanical torque gradient. Anguler 
positions are measured electrically by suitable calibration. In this applica- 
tion the microtorque, acting as a variable resistance, was mounted in the 
fuselage and attached by a linkage to the elevator (Figure 4 ). Recording 
elevator position required two leads from the model to the center of the cir 
Cle. Enameled steel flying wires were used for this purpose. From the oper- 
ator's position the leads were continued to the recording apparatus on the 
ground (Figure 6). The steady-state elevator position was recorded by 


calibrating a voltmeter to elevator position. A recording oscillograph was 
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employed to record the pulse. 





No in-flight measurements were possible in practice for two reasons. 
First, the vibration induced by the running engine was sufficient to destroy 
the electrical continuity of the micro-torque. Second, the friction intro- 
duced into the elevator system by the micro-torque linkage made control of 
the model marginal. 

Weights adjustable fore and aft were used to alter the model C.G. 
position in running the steady-state phase. Normal C.G. for the model cor- 
responded to normal C.G. for the full scale Cessna and was at 27.0% mac. 

This position was used for the dynamic phases. 

Dynamic pitch-angle data was recorded by means of a motion picture | 
camera tracking the model from the center of the circle. This camera was 
mounted on 4 heavy tripod equipped with leveling bolts and a revolving platen. 
A flash bulb in the 3rd line circuit, placed so that its flash would show on 
the film, was used to mark the instant of application of the pulse. The 
camera wes equipped with a cross hair reticule. With the camera field verti- 
cal and the camera adjusted to track in the horizontal azimuth plane, these 
cross hairs were a reference for both pitch angle measurement and vertical 
displacement. The camera speed was calibrated for a time reference during 
the transient response. 

The use of an armory was obtained for a testing site. This reduced 
the amount of stray air currents as compared with an outdoor site, and improved 


the chances of obtaining reproducible results. 
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TEST PROCEDURE 


The C.G. location of the model was found by suspending it from 
three positions having wide angular separation about the pitch axis. The 
C.G. was at the intersection of plumb lines hung from the points of sus- 
pension. 

The model's moment of inertia about the pitch axis was determined 
by swinging it as a compound pendulum. The method is described in Reference 
6 and calculations given in Appendix G. 

This C.G. location (27.8% mac) and moment of inertia (.012 slug ft") 
were used in calculating parameters for the equations of motion of the model. 

Due to restrictions imposed by conditions at the flying site (low 
ceiling lamps) and by the light weight of the model, it was necessary to fly 
it from a 37.5 foot radius instead of the 50 foot radius originally contenm- 
plated. This resulted in a higher rate of yaw, but enabled the outward yaw 
angle of the model es set by the line guide on the wing, to be reduced be- 
cause of the increase in centrifugal force with decrease in radius. The 
shorter radius was objectionable from the viewpoint of slipstream disturbances 
as the wake had less time to die out before being flown through again by the 
model. 

A. Static Phase 

No static-phase runs were made due to the above-mentioned mechanical 
difficulties. For the determination of the stick-fixed neutral point the 
model was to have been flown in steady level flight at a given C.G. position 
and velocity. The velocity would have been measured by timing a number of 
laps. The trim elevator position would have been measured with the calibrated 


voltmeter. A number of flights were to be made at this C.G. position with 
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various speeds. The speed would have been varied by a combination of engine 
adjustment and propeller change. This procedure of flying at various speeds 


was to be repeated for a number of C.G. locations. 


B. Dynamic Phase 
With the pulsing mechanism wound up, the model was flown level and 


the speed was measured. The electrical leads were connected to the opera-~ 
tor's handle from the recording apparatus, and the recording oscillograph 
started. The camera was also started, and tracked the model under manual 
control. The pulse was triggered, the pitch angle data being recorded by 
the camera and the elevator pulse by the oscillograph. The record was con- 


tinued until the stable transient died out. 


RESULTS 

No static phase data were possible because the micro-torque angle 
indicator would not record under the vibrating conditions produced by the 
engine of the model. 

A total of fifteen flights were made during which the model was 
adjusted for controllable flight and qualitative data were observed for the 
dynamic phase. These consisted of approximate measurement of speed and 
visual observation of the following: characteristics of the phugoid mode, 
response to 5, 10, and 25-degree elevator pulses, and the behavior of the 
wake made visible by exhaust smoke. 

All elevator-locked flights were above datum. Speeds observed 
were between 30-50 mph. The characteristics of the phugoid were erratic. 
Its period was observed to vary between 1.5-2.5 seconds. In some instances 


it was seen to diverge sharply. In one instance, during 40 seconds of 
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elevator-locked flight, several phugoid oscillations were seen to commence 
and damp out, under the influence of wake disturbances or other effects. 

The 25-degree elevator pulse was found necessary for a visually 
observable pulse. 

At the last pulse during ons of the runs, the energy remaining in 
the mechanism was low enough to appreciably increase the duration of the pulse. 
The result was a greatly increased transient leading to a violent phugoid 


oscillation which almost terminated the entire experimental investigation. 


DISCUSSION 
Static Phase: 

Various power settings are required during tests for elevator power 
and stick-fixed neutral point using harnessed flight. This results in changes 
in elevator power and the necesgity for evaluating thrust outputs. Thus data 
reduction is complex, and since these characteristics may be found more easily 
from wind tunnel tests it was concluded that their calculation fram harnessed 
flight tests is inadvisable. Another deterrant is external disturbance due to 


flight in the slipstream of the previous lap. 


Dynamic Phase: 

The initial intentions of this phase were to collect and analyze 
the flight data using the photographs of the time history of the pitch angle 
to plot 2 egainst time. This plot was then to be graphically analyzed to 
find the frequency response as described in Reference 11. This experimentally 
determined frequency response curve was to be compared with the frequency 
response from the analysis of the transfer function of the equations of motion 
for the conditions of harnessed flight and of free flight. Since the equations 


of motion developed in Appendix A are non-dimensional and include the harness ing 
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forces as a separate term, they apply equally well to a model or to an air- 
plane dynamically similar. Because the small size of the model prevented 
the attainment of a dynamically similar moment of inertia, the equations, as 
solved, refer not to the full scale Cessna 140 but to an airplane having all 
the Cessna longitudinal stability derivatives and characteristics (References 
7 and 8) except for the mass and pitching moment of inertia. These two char- 
acteristics, found by direct measurement of the model and calculated, are 
those of an airplane dynamically similar to the model built for this investi- 
gation. 

It was assumed from the information in the AMC report (Reference 5) 
that the stabilizing and demping effect of the harnessing line forces on the 
phugoid oscillation would eliminate this mode from any consideration in 


the analysis of flight data and would prevent it from becoming a problem in 





the testing technique. This did not prove to be the case. Although on 
flights with the elevators under control of the operator the model appeared 
stable with the C.G. at 27.5% mac; when the elevators were locked the model 
exhibited a very short period phugoid (14 to 23 seconds period) when flown 
at elevations of about 10' above datum. Furthermore this phugoid oscilla- 
tion was quite unstable, increasing to double amplitude in less than one per 
fod on numerous occasions. It was sufficiently unstable to make pulsing 

the elevator at the datum level (4%' above the floor) when flying at 40 miles 
per hour possible only at the risk of destroying the model. It was decided 
that to enable data to be taken safely a higher datum plane was necessary, 
i.e., the operator standing on an elevated platform while flying the model. 
See sketch on page 19. Since the location permitting this was not availably, 


no quantitative data were obtained at present writing. Instead, the emphasis 
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of the investigation was shifted to determination of the causes of the erron- 
eous prediction and to analyzing methods which would make the fulfillment of 


the aims of the original investigation possible. 





Sketch showing method af Yaising 


dotum plane 


In an attempt to find the causes for the discrepancy between the 
observed flights and the theory, the assumptions used in writing the equa- 
tions of motion (Appendix A) were re-examined. It was confirmed that the 
angle of the lines leaving the model in the drag direction could be neglected. 
This angle was less than 2° under the most extreme flight conditions. It 
was also found that the angle in the gravity direction due to the weight of 
the lines was also negligible. Hence the assumption that the lines extend 
in a straight line from the operator to the model is valid even in so light 
a model. This is due primarily to the short radius employed. The conditions 
that the model must be charged with 3/4 of the drag of the harnessing wires 
and that the mass of the wires and their inertia effects upon non-steady 
motions can be neglected (due to the small mass of the wire) were also veri- 


fied. Calculations supporting the preceding are shown in Appendix D. 
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However, observations showed that at the short radius flown (which was 
necessitated by indoor floor space and ceiling obstructions) and at the scale 
speed, the assumption that the circling flight would produce no aerodynamic 
effects on the longitudinal motion was not completely justified. Calcula- 
tions show that the speed ratio between the velocity at the outer wing tip 

and the inner wing tip is 1.175 for the model used in this investigation. 

For comparison, to show that the short radius used here was no more detrimental 
to this condition than the radius used at the Wright Field testing facility, 

a calculation was made based upon their B-29 model (Reference 5). This re- 
sulted in a speed ratio of 1.22 between outer and inner wing tip. Therefore 
the model used in the present investigation was less affected by this varia- 
tion of velocity with spanwise position than the large one flown at Wright 
Field on 150' lines. This speed ratio effect would produce a very complicated 
solution to the equations of motion. Therefore in the solution of the equations 
of motion it is omitted, with the reservation, however, that at least portions 
of the discrepancies between the theoretical solution for the frequency response 
and the experimentel must be chargeable to this omission. An additional aero- 
dynamic effect produced by the circling flight path is the interference effect 
of the model's own wake. Flight observations showed an intermittent disturb- 
ance which was most easily noticed by a sudden inward banking oscillation of 
the model. This was determined to occur when the model's outer wing passed 
through the slipstream remaining from a previous circuit. The banking oscilla- 
tion was produced by the fact that in counterclockwise circling flight the mt- 
board wing was being influenced by an upward component of both the slipstream 


and the inboard wing tip vortex. Confirmation of this reasoning was assisted 
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by the fact that the exhaust smoke of the engine made the slipstream wake 
visible. This wake was seen to leave the model approximately tangent to 


the flight circle with some downwash. See sketch below. 


Raqion of 
UP Word flow 





Mid e! NAG On 


circle 
Follow na cireuit 


SkeTch Showing nase effects 


Since, at scale speed the model passes the same azimuth point on the flight 
circle, and the wake appeared to widen as fast or faster than the curvature 
of the flight circle, the model crosses this turbulent area unless it is 
flown above the previous path by a reasonable distance. 

All of the foregoing has not indicated any conclusive reason for 
the unstable phugoid which was observed. It was noted, however, that in 
both the analytical solution made by Braun in Reference 5 and in the initial 
solution for the present investigation, the familiar assumption of gliding 
flight condition, namely Cp =0, was used. The solution obtained for this 
condition (Appendix C) using the equations developed for harnessed flight 


developed in Appendix A, gave as roots for the modes of motion: 
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XK, = 1.438, A = -5.35 

Nage- OL + 263; 
This indicates a very lightly damped, stable phugoid, with a period of 5.1 
seconds which damps to 1/2 amplitude in 30 cycles. This, of course, was 
completely contrary to experimental observations. It is also of interest 
to note that the short period mode has degenerated into two real roots. 
This indicates that the model was overdamped, caused by too low a moment of 
inertia. This moment of inertia discrepancy had been expected. (See Des- 
cription of Model). 

In order to evaluate the effect of thrust on the equations, the 
characteristic equations, both with and without the harnessing force, were 
re-solved using the following thrust estimate: Cp = Cp for all flight con- 
ditions close enough to the equilibrium condition such that the small angle 
approximations hold. This assumption makes Cp)*= 0 and ie O. Solution of 
the characteristic equation (Appendix C) yielded the following roots: 

d= -1.705 = -5.884 

Nay = WOB5 + .265) 
This solution indicates a trend in the right direction since it gives an 
unstable phugoid of essentially similar period to the glide condition (5.06 
seconds), but which increases to double its amplitude in 5.32 cycles. This 
corroborates the observed instability of the phugoid although it does not 
show the rapid divergence which was observed in tests. It also shows that 
the power has only a very slight effect on the short period mode. 

To check the validity of using the previous thrust assumption an 
additional solution was made of the characteristic equation (Appendix C) 


assuming that Cy was constant at the value of Ch at scale speed. This makes 
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g," = 0 and Cy > Dy a The changes produced in the coefficients of the 
characteristic equation between this and the previous assumption were so 
small that they were within the accuracy of the calculation method used to 
factor the characteristic equation. It was therefore concluded that the 
previous assumption was sufficiently precise for the size of the factors and 
for the method used in extracting the roots. 

As a further correlation of the period of the phugoid obtained 
from the equations with the period from observations the value of the period 
was corrected for flight at conditions other than at datum level. This was 
accomplished by the approximate relationship w = Cc" = Cc. 


a V2. 
Reference 1. Using the standard flight conditions and simple trigonometry, 


from 


a force of -.119 lbs per foot of height from datum was calculated, (measuring 
directions and forces in airplane stability axes direction). This is con- 
verted into a *. - correction factor by dividing by the weight and adding 
one: Ci, eC; (1 + .078h ). Therefore denoting by subscript ( ), condi- 


tions other than at datum, the following approximate relationship can be 





expressed: 
ja = = = C. eal C. a = . 
ie ai L, C.(1+.078h) } + .078h 


From the above it can be seen that Plying the model at 13 feet above datum 
approximately halves the period. Therefore the observed results agree 
reasonably with the analysis insofar as the period of the phugoid is concerned. 
Also the quality of the stability of this mode is in agreement, namely they 


are both unstable. The only discrepancy is the rapidly divergent oscillation 
which was observed but which failed to materialize from the theory. This is 


unexplained. 
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The foregoing discussion has been primarily about the phugoid 


| mode since this mode indicated the discrepancy in the initial assumptions 


upon which the equations of motion were solved. Now that reasonable agree- 


/ ment has been established between the theory and the qualitative observa- 


tions, the assumptions which produced this agreement, namely e = 0, Cy = 0 ; 
will be used in all further analysis. The phugoid mode is in itself of no 
interest in this investigation except as a complicating factor, since it is 
most markedly affected by the harnessing force of the lines (Appendix C). 
Exemination of the short period roots, on the other hand, show practically 
no difference between those of harnessed and unharnessed flight. It is this 
fact that appears to make this testing method desirable for evaluation of 
dynamic stability above the frequency of the phugoid. 

It is considered insufficient to compare the dynamics of complex 
systems on the basis of the characteristic equations alone. Such an analysis 
only considers the natural modes of motion of the system when disturbed, but 
not under the action of a forcing function. Since the normal dynamic testing 
methods involve the application of a forcing function, it is more reasonable 
to compare harnessed and free flight on the basis of the frequency response. 
This frequency response can be obtained in the conventional ee by solving 
the transformed equations of motion for the transfer function 2 as shown 
in Appendix B. The transfer functions to a step input were reduced to the 
frequency domain by the graphical logarithmic method of Reference 2, and 
illustrated in Appendix E. The frequency response, as obtained from the 
analysis of the equations of motion for both harnessed and unharnessed flight 
are plotted on Figures 7 and & . These are the curves with which the 


experimental frequency response was to have been compared. 
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Examination of the frequency response curves Figures 7 and & 
show that above a frequency of 1.0 radian per Ky, the frequency response 
for harnessed flight is identical in shape with that of the free flight, but 
the amplitude of the former is less. The fact that this difference in ampli- 
tude is a constant factor is quite evident from the logarithmic plot of 
Figure 10. On this plot there is a constant 2 decibel gain difference. 
Converting this gain difference into amplitude it was determined that the 
amplitude of the harnessed model was 80% of the correct free flight value at 
all frequencies above © = 1.0. This loss of amplitude is not unreasonable 
since the lines offer restraining forces which would tend to limit the ampli- 
tude of oscillation in response to a forcing function. 

Since the aforementioned loss in emplitude is tied up in numerous 
terms of the transfer function, it would be difficult to evaluate it directly 
from the harnessing force. However, it would appear that the change in anm- 
plitude is a function only of the mass, speed, and flying radius of the model. 
Therefore it should be possible to evaluate this percentage change for a var 
lety of permutations of these items and tabulate the result to be used in ad 
justing experimentally determined frequency response curves to free flight. 

In order to see if the transfer function of the harnessed flight 

: gave a transient response which was similar, qualitatively, to that 
observed when the model elevator was pulsed; and also to establish desirable 
frame speed for the recording camera, it was considered necessary to solve 
for the transient response. In order to simplify this calculation, a square 
pulse was assumed instead of the triangular one actually imposed on the 
model. It should be noted that a triangular pulse was chosen for the model 
to enable the gimplified method of Reference 11 to be used in converting 


from the experimental transient response back to the frequency response. 
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The mechanical details of determining the transient response from 


| the transfer function are covered in the literature. Numerous methods may 


‘be used. The one chosen for use here is the Heaviside expansion which is 


described in References 1 and 4, namely, 


Given the transfer function gi d) = F(A) 
(A) 


and roots of G() are \, teNawNae ,, none of which are zero or re- 


peated roots, 


Then = (t) = rte + Ss [coal where G(A;) = dG (A) 


This solution is the time response to a unit step input (Appendix 


| E). By using the method of Reference 4 the response to a pulse of finite 
| width can be found. This method cnnsiders the pulse as being made up of a 


| positive unit step followed at time At by a negative wit step. By super- 


posing these two inputs a unit pulse of width AX is obtained. By the prin- 
ciple of the Convolution Theorem the output is then the result of superposim 
the outputs to these step inputs. These outputs are identical in shape but 
of opposite sign and separated in time by an interval Ak. ‘These operations 
were accomplished on logarithmic coordinates and the results are plotted on 
Figures || and {2. These are time histories of the pitch angle © per 
unit elevator input. The solutions are for a pulse of .O4 seconds, .10 
seconds, and .14 seconds. They show that the pitch changes of the transient 
are of reasonable amplitude for these two pulse widths, provided that large 
elevator deflections are used. This corroborates the visually observed ef- 
fects. They also show that the pulse must be kept short or the transient P 


response due to the short period mode is negligible compared to the phugoid. 
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If a camera record at 100 frames per second were made, the pictures at each 

1/100th of a second would give ample data from which to plot an experimental 
curve similar to Figure |2.. Also, since the slower decaying of the two short 
period roots was reduced to 13% of its initial value at .7 seconds, only the 
first 7/lOths of a second would be used in the graphical analysis (Reference 


11) to obtain the frequency response of frequencies associated with the short 


period mode. 





28 


PRINCETON UNIVERSITY PAGE 
AERONAUTICAL ENGINEERING LABORATORY REPORT 181 
CONCLUSIONS . 


The use of harnessed models in circling flight for evaluating 
longitudinal stability of full-scale airplanes has been investigated quali- 
tatively. 

It was determined that this method is inadvisable for finding elev- 
ator power and stick-fixed neutral point because of the variations required 
in test power settings and because of continuous slipstream disturbance. 

The method appeared favorable for obtaining transient pitch-angle 
response to a pulse elevator input, employing a high-speed motion-picture 
camera to record transients. 

There is a definite lower limit on the size of models which can be 
used in this method, considering dynamic similarity, due to the critical 
weight requirement and the necessary mechanical installations. 

The use of an indoor flying site featuring a raised ~~ for 
the operator is recommended. This would permit initiation of transients at 
the datum plane of Ad=0 (flying wires horizontal) and would prevent destruc- 
tion of the model in the event of ensuing unstable phugoid oscillations. 

Previous to recording data the model should be flown below datum. 
To preclude effects of wake disturbance, data would have to be recorded with- 
in one lap of datum-level flight. 

If further investigation of this technique is intended with the 
Cessna configuration, it is recommended that a model of at least twice the 


size, i.e., one-fifth scale, be used. 
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APPENDIX A 


Development of general equations of motion to include Saf wbba tak forces. 
This development is essentially similar to that of Reference 5, and uses the 
same assumptions. 

Assumptions are as follows. 

(1) The lines extend in a straight line fram the operator to the 
model. 

(2) The height at which the model flies to maintain the harness-~ 
ing lines in the horizontal plane is referred to as datum height and 
the model flying at this height is flying in the datum plane. In 

this plane the wings of the model are assumed level. 
(3) ‘The bank angle of the model (A) because of the action of 
the harnessing wire guide on the inboard wing tip equals the angle 


the lines make with the horizontal. 


p ad puree — 


(4) For small displacements of the model from datum, sinad =a% 5 
cosh) = 1, and F remains constant. 


(5) The equations of motion use stability (wind) axes. 





(6) Aerodynamic effects due to curvature of the flight path 


are neglected. 
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The equations of longitudinal motion are as follows: 


=F = mV = - CokpVv’S -Wsine + Cr+ pV'S (1-a) 
=o. =-mVy = - Cu %eN'S + Weos7 + mV sin ad (1-b) 
fe 


These equations are then linearized in the usual fashion by considering 


each term: PUY am, @,---) 3 OF aN + Qt ax + dE ax + BF alad).. 
V . 


A ro. DAs 


* 
After letting Cpyp = Cy - Cr this results in 


mV +(C*pp-SV)AV a (Com%pV's ew)aa+Waoy =o (2-8) 
m Ve ~(Crpvs ~ ZrmV SY)AV - (Cia % pV*S)a% 
R 
(2-d) 


-“Wsinvy art + nav ad =o 
fz 


ie ys) - Cu, % eV Sc AV - Cma %eY Sead 


= an KpV Sead “ Crag ev Se AQ - Cm, 2 pV Scbb. =Q (2c) 
é 
In order to put the term ny in terms of the longitudinal variables, 
substitute Ab =H (where h is the vertical distance from datum), and dif- 


ferentiate with respect to time. Neglecting higher-order terms this gives: 


mV -Wsinw ot - Ci pVS4V “Cra qsaae E 4h. =o (2-b)2 


* 
AXR/dke =Nsiny where y=y, 7 AT 
However, since initial conditions specify level flight at datum level, 


m 
So =O, and am 4% = As , and from equation (l-a), -Wsinad =C64S, 


therefore equation (2-b), becomes 


mVe + ae as ax ~ CupVS AV ~C. asad + CY ar =a 


4 
1 
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Dividing through equation a and b by Mees =q4S and equation c by 


4 eV"S<=4Se and making substitutions as follows: 


mz AN ff: & pVS 
m= Vv h + Se 
ye 
v7 ps Ke hy 
me i 


| and using the operator d to indicate differentiation with respect to *4, 
| and considering ell «,¥ as being the displacements from initial conditions 


(14.e., 4%, 48 ) gives the linearized non-dimensional form of the equation: 


These are expressed in operator forn. 


(d+Co + S Cp, )m + £ Cog + tC.v =o = 


Codu + 1 Ciuda-(d>+iCod +f) = 0 a 


N]- 


mei , cl ad [da = \ cw ae 4 Cma.}d - t CLI 
h ———_ h 


+ (d* = Sade c\) = peCnyals ae 


In the analysis following, Con and C w,, will be considered small and 
‘neglected. 
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APPENDIX B 


Development of characteristic equation and transfer function S 


e 
Solution to the equations of motion (Equations 3, Appendix A) is made 


by assuming a solution of form 


AXA AX ley 


Ak 
AK = AL, O -& ro @ =a 2 


Oo =9F@ 


b 


| with the 1 Cmg de =0 for the stick-fixed condition with S- measured from 


| equilibrium elevator position. 


This results in 


Wee. d)u, + (Cox - Gigs + “ene, = © (4-a) 
C, Aw, Pee (Co + Cus 1 Cea + a, - N« C Co A+t}e, - =O (h-d) 


LCmaad * Curate, + [Omar -X je (4-c) 
TA = NW 


For these homogeneous algebraic equations to have other than a trivial 


| solution the determinants of their coefficients must vanish. 


Expension of this 3rd-order determinant gives an equation of the form 


Bees bv + bX + b.N +b, A + Ds = 0 


By defining the combination 


(5) 
Cae + Cudesm and Cited "= Cet 
A 


the coefficients of the characteristic equation above are expressed as 


? 


a 2 ~ a 
bi = Cam -Culo® Cmgy + 2 Co™Cme + £(M-Cot) + Co Cur Cua, 
= -¥ Bs a. A h 
b, = Re Cree +£(Cme + MCD”) 
uu fe hk 
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Solution for the transfer function 2 to a step input is made by 


€ 
| using the Laplace transform on the equations of motion (Equations }). 


Since all initial conditions are assumed zero for the variables used, and 
| the transform of a step function is unity in the Laplace transform as de- 


| fined in Reference 3, this results in equations 6 below. 


NIN +. (Coa - Cr} + C.8=0 (6-a) 


A 
& 


Au + |’ +(Co'+ Cua) +f la -[r CoA +f \o =O (6-b) 
[ Crnda + Cun | 4 [ Cmas \ -X]e = “LCmade (60) 


Note: Although the left sides of the equations above appear similar to 
equations 4, in equations 6, 4,4 and O are variables of ky, while in 


equations 4, “4,,0,, and 8, are constants. Solving equations 6 by determin- 


ants for 2 gives a solution of the form: 
@ 


8 (, - Gee + fol Ae | oN +O, 
Se bsrA* - Syaw + bsrX> e bin - bA > eb. 
Where the b coefficients of the denominator are those of the character- 


(7) 


istic equation and 


O, = | 


> 
a 
il 


mcac 5” 7 Cue’) 
r~ 


QA, = L\ Ce = (Com Gee CaGae || +f 


On. = Cor 
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APPENDIX C 


Since the characteristic equation is the denominator of the transfer 
function 2 (A) » all solutions will be made for the transfer function only. 
Ali solutions will be based upon the following so-called standard flight 
conditions (all angles in radians): 
Vs 58.7 ft/sec ( 40 mph) 
C.@ #232 tre WW. 
QS 


Cp2 .0616 from wind tunnel polar (Reference 8) corrected by 
addition of wire drag (Appendix D) 


e = 21.53 1b RS 37.5 ft A = 25.6 
* 1.59 rt2 Te ehh h = — ,0843 
b = 3.28 £ fe .1n19 m=  .0476 slug ft@ 
AR = 6.76 pe  .00238 Cm le -.406 from Reference 8 
C = .kg2 ft Cus 4 62 C .869 from Reference 8 
Ge ess 
Da i eTAR” 116 . 


The following various thrust conditions were considered using harnessed 


flight equations: 


Condition I -~C.~=0 constant 

% 
Condition II - C, = Cy at all « and wu (Co, = 0) | 
Condition III - Cr = Co at equilibrium o& and «, and is | 


constant at that value (Cos =€ b.): 
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These conditions produce the following changes in the components of 


Q,;, and b; : 
Condition I Condition II Condition III 
4 1.0 1.28 1.28 
Cp" 0616 0 0 
Con 16 0 6 
_ 4 .O2TT .0269 .0269 


The following terms depend upon tail efficiency ",. The Condition values 


were obtained from References 7 and 8 and corrected for h,. 


Cm. -.751 ~ 964 - 964 
ae - 149 -.190 -.190 
Cuda - 0865 - .0594 - 059% 


Solution of the transfer functions for harnessed flight gave 


Condition I: 8 (A) - _-8.91 (\'+ 2.37 N+  .252\ + a 

& > + 6.99 + G.249 rF + .TO5SGA & 5055 A + .O1 

‘ Condition II: 8 (a): -11.43 (A+ 2.31LA + 1h 

6 (A“ + 7.577 A> + 10.105 A 7 W125 \ +. 702) 
Condition III: 9(,)- _-12.43 (A+ 2.32 + «142 

6 (NM + T-5TTA® + 10.152 AY + .3023A + JOC) 

Factored forms of the preceding are: 

Condition I: 9(,)\~ _-8.91 (A+ 20262) ( + 084) (A+ 2,256) 

6 he O327)(A + LA3E)(K + 5.53) (A + 002A + .2027°) 
Condition II: 8()- _-11.43 (A+ 065} ( + 2,2h5) 

| & N+ 1.705 N+ pie + .O12 + .2 


Condition III: = (A) = Same as Condition II 
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Condition II was established as the most proper condition and the trans=- 


| fer function was re-solved for unharnessed flight by computing the coefficiemts 


'for f=0. 


| Condition II+- 9 an A) = -11 .43 (A a en te +_ 2.298 
| Free flight: N+ lll N+ 507d) (A* - .OOLGA + 1.113 











APPENDIX D 


| 


. 






Calculation of drag effects upon the harnessing wires. This solution 


ssumes no outward yaw of the model. Axes are as shown. 


| 





Vv* 
T e# line tension ® F = 7 


Cp ®= 1.2 two-dimensional flow Cp of cylinder at low Ny based 
on frontal area. 


Sw °° area per foot of wire normal to velocity 


= .0007 for .008" wire 


-OOl for .012" wire 


q, = 1/2p-5° where ~~ is velocity of wire at station x 








Taking the summation of forces in the drag direction, assuming T remains 


constant over the element, 
ZFY = Tsin®O@ - Kxt - Tsin(O + a2 | —®) 
Ox 


Tswoe - Kx* - TsinwOcos dO -~ Veos 8 Stn dO = 6 


Ax dx 
cosO=Z!l ,cosdea =!1, Sindd =de = cl‘y 
d x K dx Ae 
dA x* 
Integrating gives 
dy = - K x? * Cc 
ax aT 


Integrating again gives 


pees Wee yw [> 


ut 1. 2a 
Solving by constants for end conditions: 
D=0 
3 
= 2 
Therefore (ae 
3 
mo 8 Eliz «CK R 
-.. 2 
S) 7 
i = -30e\ | 
MmaCof/V\c . 3S 
K =p S°(X)S = 38,8 
eh, «8 ee eee = - 52 S 
Sm 


For tension shared by all wires equally (most conservative condition) 
Ong = (-23.72)(.00234) 57.3 = 1.840 
Drag of wires: 


ad 2.5° =: 350 ee! 
Operater ~~ 
Control handle 





D = ADdx <= K x*dx 


UO 
£ 
’ 
sh 
S 
= 
c) 
ee 
Ls 
oe 
r 
WY 
| 
SR 
om) 
\O 
E 
rH 
O 
ty 
o 
S 
= 
4 
@ 


ae es 
=(4)Do, = 1244 2. 
Determination of proportion of wire drag adding to airplane drag: 
Centroid of wire drag: 





Xx ' i+ QR. 
==— ae tlt 3 
D D Veins = Bf, Kx Ax 
31.5 uy. 37-5 
= 3 Oe eo 
D NMxdx a qj Bile 
2.8 
3 32.5 
bur DO= Kx \ 
3 32:5 
Therefore \< a 31.¢ 
Ko = q : 2 7 v1.8 
i< x3 as : 
a a1.S > 





Therefore the wire drag centroid is located 3/4 of the distance along 
the radius from the operator. 
Dwires © 3/4 (.144) = .108# 
Converting to conventional coefficient: 


Co,, ® 1068 = .0166 


6.52 
Total model drag coefficient: 
Cy F Te © ee us © vuine 
Tetal drag: 
D = Cpo aS = -O616(6.52) = dor 
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er 


Calculation of Line Tension: 


Centrifugal Line Force T, = “VY _ WY 
> 


Side Force @ S° of rc -<Cy>Cu,gG = .008x5 = .OhO 
Y. = CuoqgS = (.040)6.52 = .26 lbs 


Total Mine Tension = 4.64 ibs. 
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APPENDIX E 


Frequency response calculation of model in harnessed flight, C+ = Co: 


-rst (t+ ty ee) | 
Meee eNO ee 


| 35 . 
(i+ Seat + ee) LCRa) - 18% + 4 


Components of the Lm-Angq plot are given Figure 4 and are added to 


6 ee = 
ra ) 





obtain the complete frequency response in Figure !0. This result is re- 


plotted on cartesian coordinates in Figures 7 and 8 . 


ee = 2.31 
2.6 Lou... 2.31 2=1.S8 db 


lus 2 (ys) = 2OLon,. 24.3) + 2.0 Ley.e J +(32) + no Leyes f) +( Ae Y 
- 20Lkoyo fi 4/32) | 20feyn fy 4/32 ~ 208. fie (d2 ) Paci 


d of 2 (a0) = Low’ Bw + Weave WwW een Tae 
O65 2Z.2zads S.e84 
-Se > _ kan (-.1etw) 
{.7eag a ee ky 


rr) 








; =) 
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APPENDIX F 


Transient response by Heaviside method: 
Cw * 
2 (x) = aS) Ss ays aes 


QCo) “i Sy IAS) 
Harnessed flight, C,=Co : 


" “11.430 + 2.31\ + .1459) 
G&A & TTT AS + LOLLOD AY + ALAS + «702) 


Plo) _ -12.43(.1459) = -2.376 
Qo) (02 


A, = -1.075: 


P(A.) = 2.907 - 3.939 + .1459 
= . 896 


Q(A\= ¥AF + 3(7-577) A. © 2(10.165)A, + 14125 


= -19.82 + 66.079 - 34.66 + 4125 
= (2.01 


OCA) _ _g86(12..43) @  -. 7044 
AN, Q'CA,) -1.075(12.01 


rx = ~5.,885: 


P(A.) = 34.62 = 13.592 + » 1459 





= 21.174 
Q(A.) = -814.8 + 786.95 ~ 119.96 + .4125 
= -1h7.40 
KPC) — -11.43(21.174 
AQA) =e any. Ey * 
2 (x) = 2,376 - Token ts ee .2790 €?°  BBh 


(Short. period mode only) 
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Long period mode of harnessed flight, ee one 
n° B(e\= 4.1297 Asin (ll. + 135:5°) 


The complete solution for this condition is plotted in Figure !2. 






Short period transient, free flight, Cx =n: 
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Sample calculation of Moment of Inertia Ly 


Rass eee” .(e + Ve)d* Te 


oe 
Model weight =~ 1.53 lb. 
Pendulum length <= 1.23 ft. 
Swinging gear weight = .0331 lb. 


Swinging gear pendulum length = 1.053 ft. 





Q 

V= Measured period of oscillation = 1.335 sec. 
N= Model volume = .216 ft. 

P= Air density = .00238 slugs/ft.- 

= 


= Gravitational acceleration = 32.2 ft _ 
sec. 


I. = Inertia of swinging gear 
Pees 2 
Yn” 


where { = measured period of oscillation of the gear alone = 


000975 





Substituting in the above equation for ly, : 








a 
- if Ee 
7 rok * - 
fn 
ft 
7 






i 6 6} 
ase° 
508° 
h . 
y Y W" 38° gyatic 
ad = . - — om I ee 
AREAS — 
WING (INC FUSFLAK) 19929 SOFT 
ST.MtLiZER 1@ 6eq 
Fim 6668 
ELEVATOR (PLUS TAB) 9 @B = 
RUDOER 4.752. « 
SILEROA (TOTAL) 14 08 ve 
FLAPS a7" 
SLEVAYOR Tae ry tr 
ANGLES — 
INCIDENCE STABILIZER -Ef° 
INCIDENCE wine ro the 
OME ORAL ree 
) AUDOER TRAVEL ¢1@ °° 
~ , &@watoa Teaver 220 3! 
ewe % ss 4 
AILERON VAaveL 4 Re 
de e 
. ERMME CONTINENTAL Bs 


“a TT) 
are. 2075 
" FROP. 818 70° (max.) 
‘ \ 
: WEIOHTS 
cupry ere — € REBAR SPAR 
ernose 150 = a 


weeene 


Al a 
~<a 36 


¢ From’ Par 










Teri ws ae oe ff) UE oe & a8 Gs 128 
ro» ome an 








ae Gmhe 2. 2 Owe Geaw w *¢ aves OF) USE 
= LL. @ Gam - aD @. 6 Saeeme Ores 


NEUTS 


toe oe 


wiv Ww 





> 
3 
ta t 
i> = = mM 
a - —rtreomir 
b ae: ts Frits | 
sh gd “meemelse 
Bitte! rare | f t 
t Py 32h iz078 . : 
Bee 3 fh om Pe j : 
Eiki agpects | 
zs o° oy ‘ a 3: j 
Vice 3 “ $3 ~ 
—_— pes . . : 
(2 es | 12 ene Spee 


ey Se eee 


SE a Ca a PS RIM SPP <P = — 


ries. 14& 2 





FIGURE 1. - DISASSEMBLED MODEL 
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FIGURE 2, - COMPLETE MODEL 
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FTGURE 4a. - SHOWING PRIMARY CONTROL WIR&S AND COMTRCL HANDLE 
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DETAILS OF MEASURING 









APPARATUS 
Micre- Trager wire and Sead baite cy Vor cane 
Te vgue £ oo Aq wicks Nitin oe al tap Yetistor 
: a ] |} 
F a pa 22.5V 
Flash bulb eet 
«Cam 
S : | 
AN? : 
> 240 Su : : 
Figure 5 7 
Circuit Biuaram RECO ding 
ee ee Oscitrs araph 


rig uyre b 
General Ar rangement 


/ Flying WITES 





Teraqer wWwive ~™ 


oY ae oN ~ \n 

iif 1 \iyo.cu \ 

ei | ee ‘a A Taeecu 
Seno — oe. _ 


ee Aaol Fie? 
+ 


= Ss 


s 





Figore ri 


f . . - | | \ | : 4 

fe oe | “i | rey rq 

oS Pes ed eo 
- 3 i } . a - 


(bd ect 2h a ecw ee Cry —s 
| . ; mi 


= 
rr 
ve) 





ee 








: 
— , ey 
—|--+4 4 dong tu 
| 4 xIpus ayt40| a uy - Wy no 
ae | = | 
| 2 BE EME oF oe 
a Yap ee ee 
| | 
at an Tor a eid ae 


ase h 


a 
Lian maT tis nella way 


| 
| ; 

| 
: | | : e 
z ! Pee ae | ee . | ‘e) 
| ~ ap Suse - ae) oo a a ae aA Cl CeCe ieee Fe ei bik) Me we = 

j i | ) : 

| | 


a 
~-~;- 
| : 
~~, 
| 
ma 
| 
ene VK = ae 
i 
' 
-— 
' _ 
_- 
I 
{ 
| 
r 
| 
a 
t 

















Figure & 


Ji 

+t 

_ = ‘cee 
— ; { 

{ | 


addy ‘god Sug 


. 
4 
5 
ed 
e i 
‘ 


- ee 


ors | | _ oe . | 
sop ataypiis eoba ete ope 
= ee ee mh a ee a st a eee 


Ty 
ay ' 


1 $ 
a er er ~a 
? ‘ 


4 
- ry = 





— ; ae eT a =" Semen ee 


—epaentts  le 





-4- 


Vis 
TERS 


oe 

ci Cc 
Mg 

ma be doo 


ale 
aS 


+t 1 
f 
a7 


es 
' IMA OD 
| 7 ss 
| ret Po 


| | 
: doer 
: 
| 7 


oh — 
4 we 


ee. 





~~ 


ae 
LC 


— = 
—4— 














Figure it 


ort fm Oo MK~ YO} WH z - ~N ~n om 1} WwW =z er) ~N << 


i 


LF 
ifn 
oe i i) 
Hi 
if 
Lid 
sats He 
| 
“4 
hig 
ills 
Hl 






a Hee 


if 
siiltilt poegaestai suite 
mone | ey 


_ 


pitas 
deed tees 
g EGSce fees iste fess 


i 
Beer 


cea 


ih 


ae 
Ha 
EE gE 
htt 


BSaLE 


oes SESSee tii tiie 


tenge Hee ewer Be e 


aging es 
HU ant sie ‘ate ae 


BAe I a ul pia ttroenes 


eees4 gees wets 
i 


Se 


nes one a a 


e aeEes 
ee 









6 


9 
Sead posed sn HG 


ttt 
ines HHH er 
Hd Ure astaa iil iit 


ie 
ee 
ae "iiss 
gees 
| Ane 


od af Gna ees bss snag seme diate SHEE ee SESS = phiiette +P 
yey EESSa Besa Carte: SEREEEES jespERESS SG Hct - 
EES ils! 7aa) nN Eaas SESE SESS itt ee iit 7 
i i SSS sh Sse Bless (lec ss SSS SSS syecses= 
fant Bot enemas SSesseomeue suietie ae 
wae Ss 
init 


aa 
4 







tem ses srs CRaSE cattarteoaceaun °o eae ates ae il rae tT 
0 cao a, Senta Semen ca gE 
a aS Eee 
2 apg Stee ig Un UNS conse Tes 
pi ] + | @ sti t So Pe ae stan, 

Be 5 NE Ht ais 

6) . 
th 


7 






7 ins 
TEU gd ESET Egat 
et 


St sen 
reas 
ou feels eerrppee | Pd 





5 
deg Sees : 
Ein lid Old Sr bse Sees BEES EEET ec 





o7 
ety? | Poe 






eeguaes | 


F qure \2 





rs: ae le Te: 
1} Leal ; 
et bpd 
Fag Gal — Tk ; 
Jeo : 


35:5at)| 





SA 
SI HEE “8. 


on 
a 


AS 
bo fb 








NS€ 
t 





| 
TT ae 
a . = J . . 
ab , 
' bose 
> 
ve 1) 








ont inuk 


hc 
le 





2k 
al 
i 


202-6 927 
mh 





{ 

1 
cent aa ease 
SE 


Respo 





t 
: 
yk 


| 





sett Flial Tf, Cr. 


iz 


: 


me 


co 


Sat sint 
| 
_ 
[a 
| 





a . 
{1} — 
[am 





War 


osc\i 


a4 SinB5.5 
a 
nds 
— 


—+2 


- Seto 


3 


ym o 


> 





eval F 
ze/b2itig 











a 





i ia ah -— FD =i 
2 > . 4 
: - t> Sad 5 _ “+ oe 4 ‘ 
ic ae z as ie pe ; H 
a, ee ae a ba Poi- ta : ae 
: = . < ° a oa a ET -+_ ——= -—| = — ual | ed 
‘ 1 « / © oa t 





Bear Uae ) 





+ 


. 


ey 














Thesis 16269 
M93 Munk 
Irvestigation of the feasibility 
of using models in harnessed circling 
flight for the evaluation of air- 
plane longitudinal stability 


J ar 
7. al Tosteradnaty S 
yer, Ur tye Rad Ue 





